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A system's temperature can be denned in terms of its constituents' instantaneous positions rather 
than their momenta. Such configurational temperature definitions offer substantial benefits for ex- 
perimental studies of soft condensed matter systems, most notably their applicability to overdamped 
systems whose instantaneous momenta may not be accessible. We demonstrate that the configura- 
tional temperature formalism can be derived from the classical hypervirial theorem, and introduce 
a hierarchy of hyperconfigurational temperature definitions, which are particularly well suited for 
experimental studies. We then use these analytical tools to probe the electrostatic interactions in 
monolayers of charge-stabilized colloidal spheres confined by parallel glass surfaces. The configu- 
rational and hyperconfigurational temperatures, together with a novel thermodynamic sum rule, 
provide previously lacking self-consistency tests for interaction measurements based on digital video 
microscopy, and thereby cast new light on controversial reports of confinement-induced like-charge 
attractions. We further introduce a new method for measuring the pair potential directly that uses 
consistency of the configurational and hyperconfigurational temperatures as a set of constraints for 
a model- free search. 



GENERALIZED TEMPERATURE 
DEFINITIONS 



A variety of thermodynamic temperature definitions 
complementary to the classic kinetic definition have been 
derived recentlyfi], 0, [j| . The most general form, proved 
in Ref. 01 and 0, is: 
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where T = {T lr -- ,T 6N } = {qi,--- ,q 3N ,pi,--- ,p 3N } 
is the set of 37V generalized coordinates qi and 37V con- 
jugate momenta Pi describing the 6iV-dimensional phase 
space of an TV-particle system in equilibrium, and angle 
brackets (• • •) represent an ensemble average. The sys- 
tem's Hamiltonian, H(T) = K({pi}) + V({qi}), consists 

of the potential energy K({pi}) — Y^i=iPi I '(^ m ) an d 
the conservative A^-body potential V({qi}). The most 
noteworthy feature of Eq. is that B(T) can be any 
continuous and differentiable vector in phase space 0. 
Choosing B(r) = {0, • • • , 0,pi, ■ ■ ■ ,Psn} yields the fa- 
miliar equipartition theorem, 
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where rrij is the mass of the j-th particle. Choosing in- 
stead B(r) = -W({g;}) yields 
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which depends only on the objects' positions and not 
their momenta. This has come to be called the configu- 
rational temperature. 



The configurational temperature's independence of the 
particles' momenta has important ramifications for ex- 
perimental studies of systems such as colloidal suspen- 
sions whose configurations are easily measured but whose 
momenta are not. The experiments described in this Ar- 
ticle exploit properties of the configurational tempera- 
ture to obtain new insights into the interactions between 
charge stabilized colloidal particles. When combined 
with thermodynamic sum rules, this formalism provides 
previously lacking thermodynamic self-consistency tests 
for measurements of the particles' effective pair poten- 
tials. The same formalism also can be used to measure 
pair potentials in soft-matter systems directly, thereby 
bypassing questions of interpretation raised in previous 
studies, and yielding comparably accurate results with 
substantially less data. 

Section[n]provides an overview of several consequences 
of Eq. These are extended in Sec. lIIII to a hierarchy of 
hyperconfigurational temperatures that lend themselves 
particularly nicely to experimental studies. Section IIVI 
provides an alternate foundation for the entire configura- 
tional temperature formalism in the classical hypervirial 
theorem. Consistency among the myriad temperature 
definitions is possible only if the assumptions underlying 
their derivations all are satisfied simultaneously. Apply- 
ing these definitions to experimental data, as described in 
Sec.0 therefore probes the nature of the system's inter- 
particle interactions. These data were obtained from 
digital video microscopy measurements on monolayers 
of charge-stabilized colloidal spheres dispersed in water 
between parallel glass surfaces. Not only does this sys- 
tem lend itself naturally to computing the configurational 
temperature, but the results also help to resolve a long- 
standing controversy regarding the nature of charged col- 
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loids' interactions in confined geometries. These results 
are summarized in Sec. IV1I 



II. CONFIGURATIONAL TEMPERATURES 

The most general configurational temperature defini- 
tion, Eq. ©, requires knowledge of the full TV-body po- 
tential, which typically is not known for experimental 
systems. A more experimentally accessible form emerges 
if the particles interact via pairwise additive potentials, 
Uij{vi — Yj). In this case, we can interpret gradients of 



the iV-body potential, 



W({fc}) 



d 

dqi ' 
(Fi,. 



d 



' dq 3 N 



JV-l N 

(4) 



as components of the net force acting on each of the parti- 
cles due to their interactions with their neighbors, where 
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is the force acting on the 7-th particle. The temperature 
then may be written as 3] 



iF = 
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It makes sense that the temperature should be reflected 
in the instantaneous distribution of forces because ob- 
jects explore more of their potential energy landscape as 
the temperature increases. 

All definitions of the temperature derived from Eq. 
involve approximations of O (1/N), and so only are valid 
in the thermodynamic limit, N — > oo. Dropping addi- 
tional terms of O (1/-/V) in the derivation of Eq. © can 
be justified for systems with short-ranged interactions 
and leads to another equivalent temperature definition 
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Neglecting further higher-order terms yields still another 
form [5j, 
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virial temperature for s = 1. By the same token, we pro- 
pose that B(T) = {Ff, ■■■ , Fg N , 0, • • • , 0} with s > 0, 
yields the set of "hyperconfigurational temperatures" [j| , 



k B T h 



of which is equivalent to the standard configura- 
tional temperature. Here, Fi is the magnitude of the z-th 
element in the set of 3iV components of the forces on 

(s) 

the N particles. Because Fi is non-negative, T fl is well 
defined for any positive real value of s. Negative values 
of s would yield diverging temperatures because at least 
some of the Fi will be vanishingly small for any system 
substantially larger than the range of interactions. 

A simple example motivates introducing this new hier- 
archy of definitions. If, for example, a system is charac- 
terized by Coulomb pair interactions, u(r) = 1/r in d = 3 
dimensions, each term of the denominator, V^u(r) = 
r l ~ d d r (r d ~ 1 d r u(r)), of Eqs. © and JJJ vanishes. Conse- 
quently, the associated configurational temperature def- 
initions in Eqs. ©, J7J) and (jSJ diverge unphysically. 
The hyperconfigurational temperatures, by contrast, are 
still well defined with d x F x = (1 — 3x 2 /r 2 )/r 3 , and 
F^dxFz + F^dyFy + FI- x d z F z ± for s + 1. Con- 
sequently, the hyperconfigurational temperatures should 
apply to any system whose pair potential is continuous 
and differentiable. This suggests that they will be useful 
for studying systems whose interactions are not known a 
priori. 

Additional useful results emerge for systems such as 
colloidal monolayers whose interactions are isotropic. In 
this case, the Cartesian coordinates may be analyzed in- 
dependently 
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Setting B(T) = d x V{{ qi \) or d v V{{ qi \) in Eq. JTJ leads 
to same results. We will refer to the two terms in 
Eq. HlHjl. as well as analogous results for other tempera- 
ture definitions, as the Cartesian components of the con- 
figurational temperature, T x and T y , respectively. 



III. HYPERCONFIGURATIONAL 
TEMPERATURES 

Choosing B(r) = {<?!,••■ ,q^ N ,0,--- ,0} with s = 
1, 2, 3 ■ • ■ in Eq. (JTJ yields a hierarchy of so-called hyper- 
virial temperatures [6|, which reduce to the Clausius's 



IV. DERIVATION FROM THE HYPERVIRIAL 
THEOREM 

A remarkable consequence of Eq. JQ) is that any vector 
field B(r) that depends only on the 3iV configurational 
degrees of freedom gives rise to a functionally distinct but 
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thermodynamically equivalent definition of the tempera- 
ture that depends only on configurational coordinates. 
Here we show that this insight emerges transparently 
from the hypervirial theorem 0, 0, and that Eq. Q 
itself can be derived from this starting point. 

The Hamiltonian equation of motion for an arbitrary 
dynamical variable, f(T,t), is 
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dt 
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is the Poisson bracket of / and H(T). Provided that 
(1) / does not depend explicitly on time and (2) f(T) 
remains bounded in the course of time, the time average 
of Eq. (|llf> yields the classical hypervirial theorem, 
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If we restrict / to be a homogeneous first-degree func- 
tion of the momenta, 



f(T)= Pj Q({ qi }), 



(15) 



while still allowing Q to be an arbitrary function of co- 
ordinates, the right hand side of Eq. I|14|) becomes 
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The coordinates and momenta are statistically uncorre- 
lated in equilibrium, so that 
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where terms with i j vanish. Now, (pi^pj = k B T 
is the standard kinematic definition of the temperature. 
Substituting this, Eq. JTBJ) and Eq. JTHJ into Eq. l[T3)l 
yields 



k B T = 



dV({ qi } )\/dQ 
dq } 



(19) 



which is a special case of Eq. |QJ when B(T) = B({q , j}). 

Hirschfelder 6] chose Q — qf to obtain a hierarchy 
of hypervirial temperatures characterized by the index s. 



Choosing instead Q = Ff yields the hyperconfigurational 
temperatures in Eq. ||5J. 

A similar line of reasoning provides a straightforward 
derivation of the most general result, Eq. QJ. We be- 
gin by choosing Q = qi in Eq. (|19|) to obtain k B T = 

Qig^^, which is equivalent to the Clausius virial the- 
orem, ksT = (rjFi), for systems with pairwise additive 
interactions. Next, we choose f(T) = qiP({pi}), where 
P is an arbitrary function of the momenta, and substi- 
tute into Eq. (|14(l . The following steps are analogous to 
those used in deriving Eq. I|19|l and yield 
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Combining this with Eq. (|19|1 yields 
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where, once again, we have exploited the statistical inde- 
pendence of the coordinates and momenta. Because this 
holds for any choice of Q and P, it holds for any sum of 
products of the form ^ m Q m P m . Thus, Eq. (IL'l'li is equiv- 
alent to Eq. Q for any choice of B(T) whose components 
can be expressed as Taylor series in the coordinates and 
momenta. 



V. APPLICATION TO COLLOIDAL 
DISPERSIONS 

Because the positions of microscopic particles, such as 
atoms, are hard to measure, experimental applications 
of the configurational temperature have been relatively 
limited |5( . Colloidal dispersions differ from many exper- 
imental systems in that their constituents' motions can 
be tracked quite easily. Unlike most macroscopic model 
systems, they can be thoroughly equilibrated through 
the particles' intimate contact with the suspending fluid. 
Furthermore, an effective potential between two parti- 
cles often can be defined. For this reason, colloidal 
dispersions offer an almost ideal experimental test-bed 
for studying the configurational temperature and related 
ideas. These ideas, in turn, provide valuable new tools 
for assessing the microscopic state of this important state 
of matter. 

This section describes methods for estimating the con- 
figurational and hyperconfigurational temperatures of 
colloidal dispersions from digital video microscopy mea- 
surements of their microscopic dynamics. Emphasis is 
placed on how to account for inevitable experimental er- 
rors in systems of limited size whose interactions may 
not be fully characterized a priori. In particular, we an- 
alyze the structure and dynamics of dilute monolayers of 
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charge-stablized colloidal spheres confined by planar sur- 
faces, on the one hand using convergence of different tem- 
perature definitions as a thermodynamic self-consistency 
test for independently measured effective pair potentials, 
and on the other as a new approach to measuring the 
pair potential directly. 



A. Interactions in charge-stabilized colloid 

Considerable attention has been focused in recent years 
on colloidal interactions, particularly in light of exper- 
imental observations that challenge theoretical predic- 
tions. For example, charged colloidal spheres of diam- 
eter a and charge number Z dispersed in an electrolyte 
interact with each other directly through their Coulomb 
repulsion and also indirectly through their influence on 
the surrounding distribution of simple, atomic-scale ions. 
Poisson-Boltzmann mean-field theory predicts an overall 
screened- Coulomb repulsion of the form 

0u{r y =Z ^ B ^^^t^l, (23 ) 
(1 + if) 

where = fc^T is the thermal energy scale at absolute 
temperature T, r is the spheres' center-to-center sepa- 
ration, and Xb = e 2 / (AirelcsT) is the Bjerrum length 
in a medium of dielectric constant e. If the electrolyte 
has a total concentration c of monovalent ions, then the 
Debye-Hiickel screening length, kT 1 — l/^inXsc, sets 
the range of the effective electrostatic interaction in the 
mean-field approximation. This is an effective interac- 
tion because it results from an average over the simple 
ions' degrees of freedom. When viewed in this light, it is 
not surprising that measurements might differ from pre- 
dictions based on Eq. (|23|l . More surprising is that like- 
charged colloidal spheres appear to attract each other 
under some circumstances, in qualitative disagreement 
with Eq. J221- One of the goals of the present study is to 
apply the configurational temperature formalism to re- 
solve some of the outstanding questions regarding these 
anomalous like-charge attractions. 

Unlike atoms, which travel ballistically within the 
potential energy landscape established by inter-atomic 
interactions, colloidal particles are immersed in a vis- 
cous fluid that randomizes their trajectories over in- 
tervals longer than r = [3mD, the momentum relax- 
ation time. Given a typical colloidal diffusion coefficient 
D m 1 ^im 2 /sec, r » 1 psec. Consequently colloidal sus- 
pensions' microscopic temperatures are not easily mon- 
itored with the usual kinetic definition of the tempera- 
ture, Eq. The configurational temperature, by con- 
trast, can be measured from snapshots and so provides 
an ideal alternative. 

The fluid also acts as an intimately coupled heat bath 
whose heat capacity vastly exceeds the colloidal par- 
ticles'. Consequently, the dispersion's thermodynamic 
temperature is all but guaranteed to be the fluid's, which 



is readily monitored with standard techniques. It is natu- 
ral, therefore, to compare estimates of the configurational 
temperature based on microscopic dynamical measure- 
ments with this bulk thermodynamic temperature. 

Our samples consist of negatively charged silica spheres 
a = 1.58 pm in diameter (Duke Scientific Lot. 24169) 
dispersed in water and confined within a slit pore of 
height H formed between a glass microscope slide and 
a cover slip. The glass surfaces also develop large neg- 
ative charge densities in contact with water 0], which 
repel the spheres and prevent them from sticking under 
the influence of van der Waals attraction. Silica spheres 
are roughly twice as dense as water and sediment onto the 
lower wall in a matter of seconds. The low-concentration 
samples used in this study thus form a dilute monolayer 
once they reach equilibrium. Reservoirs of mixed-bed ion 
exchange resin help to maintain a total ionic strength 
around c = 5 x 10~ 6 M in the 1x4 cm 2 visible sample 
area. 

The hermetically sealed sample is allowed to equili- 
brate at ambient temperatures on the stage of a micro- 
scope. The particles' motions are imaged with a CCD 
(charge-coupled device) camera and video taped at 30 
frames/sec before being digitized. Standard methods of 
digital video analysis [T(| identify the particles in each 
video frame and report their locations in the plane with 
a resolution of 30 nm. The resulting distribution, 

N(t) 

p(r,t) = J2S(r-r j (t)), (24) 
i=i 

of N(t) particles in the field of view at time t provides de- 
tailed information regarding the particles' dynamics un- 
der the combined influences of random thermal forces and 
their mutual interactions. Distilling this information into 
an easily interpreted form requires further analysis. 

One of the most commonly used tools for analyzing col- 
loidal microscopy data is the radial distribution function, 
g(r), which is computed as 

g{r) = \ (p(r-r',t)p(r',t)), (25) 

where n = (p) = N/A is the areal density of particles in 
a field of view containing N — (N(t)) particles. Angle 
brackets in Eq. I|25|l denote averages over both angle and 
time. A typical example appears in the inset to Fig. ^ 

B. Colloidal interaction measurements 

For a system with radially symmetric pairwise-additive 
interactions, the radial distribution function can be re- 
lated to the pair potential, u(r). As a starting point, we 
introduce the potential of mean force, 

w(r) =-k B Tlng(r), (26) 

which reduces to the pair potential in the dilute limit: 
lim„^ou>(r) = u(r). Crowding at higher concentrations 
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induces many-body correlations that appear as oscilla- 
tions both in g(r) and also in w(r). Such oscillations 
arise even in systems with monotonically repulsive in- 
teractions, and could be mistaken for structure in u(r). 
Reducing the concentration to avoid many-body correla- 
tions often is not practical given the competing require- 
ments of tracking particles accurately (which favors high 
magnifications [ill) and amassing adequate statistics in 
g(r) (which favors a wide field of view and thus low mag- 
nification [l2|]h 

No exact closed- form relationship between g(r) and 
u(r) is known for systems at finite concentrations. Fortu- 
nately, Henderson's uniqueness theorem guarantees 
that any trial potential that reproduces g(r) uniquely 
describes the system's interactions, provided these are 
indeed pairwise additive. In light of this, two comple- 
mentary methods for inverting g(r) to obtain u(r) have 
been introduced. The inverse Brownian dynamics [T^j 
and inverse Monte Carlo techniques compare the ex- 
perimental radial distribution function with results ob- 
tained from numerical simulations based on experimen- 
tal conditions and trial potentials until a match is ob- 
tained. The liquid structure inversion technique 16] uses 
the Ornstein-Zernike (OZ) integral equation |l7| to cor- 
rect for many-body correlations directly. This is the ap- 
proach we will adopt. 

The OZ equation describes the evolution of many-body 
correlations from a hierarchy of pairwise interactions. 
Truncating the hierarchy yields approximate but ana- 
lytically tractable relationships between g(r) and u(r). 
Among these, the hypernetted chain (HNC) approxima- 
tion is found to be accurate for "soft" potentials while the 
Percus-Yevick (PY) approximation is more accurate for 
short-ranged interactions. For two-dimensional systems, 
the pair potential can be evaluated in these approxima- 
tions as UM 



u(r) = w(r) 



k B T n I(r) 



(HNC) 



k B T ln(l + nl(r)) (PY) 
where the convolution integral, 



(27) 



7(r) 



) - 1 - nl{r)} 



■r|)-l]dV, (28) 



is solved iteratively, starting with I(r) = 0. 

Provided the particles' concentration is not too high, 
and that their interactions are pairwise additive and 
neither too long-ranged nor too short, and that g{r) 
is known with sufficient accuracy, both approximations 
should yield equivalent results for u(r). Unfortunately, 
consistency does not guarantee accuracy |12J, and addi- 
tional checks are necessary to ensure that pair potentials 
obtained with Eqs. H27fl and H28fl are meaningful. This 
is all the more important if such measurements as the 
example in Fig. ^ run counter to long-established theory 
|l9j| . In this particular case, the 0.3 k B T deep minimum 
at v = 1.5 (7 is inconsistent with the purely repulsive 
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FIG. 1: Inset: Radial distribution function g(r) for a mono- 
layer of a = 1.58 /im diameter silica spheres sedimented onto 
a layer at areal density na 2 = 0.0684 above a glass surface in 
a parallel- plate slit pore of height H = 9 /im. Main plot: Ef- 
fective pair potential obtained from g(r) using Eqs. and 
l)28|l . The smooth dashed curve is a fifth-order polynomial fit. 



interaction described by Eq. (|23|) . Such anomalous in- 
teractions have been reported before in confined colloidal 
monoalyers [Til Il6l IT9I EoT ] , but their origin has remained 
unresolved for more than a decade. Is this minimum a 
real albeit unexplained feature of the charged colloids' 
pair interaction, or is it an artifact of the measurement 
technique? 



C. Temperature-based consistency tests for u(r) 

Because hyperconfigurational temperatures depend 
sensitively on particles' interactions, they constitute a hi- 
erarchy of thermodynamic self-consistency tests for u(r). 
The first successful application of the configurational 
temperature formalism to colloidal interaction measure- 
ments was reported in Ref. p|. Here, we provide a more 
detailed description of these tests' implementation and 
compare their results with those from thermodynamic 
sum rules. Together, this suite of independent tests con- 
firms that pair potentials extracted from digital video 
analysis of confined charge-stabilized colloidal monolay- 
ers can be both thermodynamically self-consistent and 
accurate. 

The flow chart in Fig. [^outlines our strategy for check- 
ing and optimizing pair potentials. Given a candidate 
pair potential, we calculate a variety of configurational 
temperatures and compare these with the known bulk 
temperature at which the experiment was performed. 
Equation l|27|) conveniently yields u(r) in units of k B T, so 
that configurational temperatures derived from u(r) are 
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FIG. 2: Flow chart for applying configurational temperatures 
to colloidal interaction measurements. From the measured 
positions of particles, p(r, i), we can calculate the radial distri- 
bution function g(r) and, from this, the candidate pair poten- 
tial, u(r). p(r,t) and u(r) are used to calculate the configura- 
tional temperatures, while g(r) and u(r) are used to compute 
the sum rule in Eq. 132H . Thermodynamic self-consistency in 
both tests suggests that w(r) accurately describes equilibrium 
pair-wise interactions in the system. Data were obtained for 
a sample at H = 18 /im. 



automatically normalized by the bulk thermodynamic 
temperature, T. The condition for thermodynamic con- 
sistency therefore reduces to T con fi g = 1 in these units. 
Successful collapse of the hierarchy of configurational 
temperatures to the thermodynamic temperature sug- 
gests that all of the conditions assumed in the configura- 
tional temperatures' derivations have been met and that 
u{r) accurately describes pairwise additive contributions 
to the free energy of a homogeneous isotropic system in 
equilibrium. Conversely, a failure to converge suggests 
that one or more of these conditions has not been met. 

Identifying departures from equilibrium is particularly 
important in experimental systems where subtle fluid 
flows or environmental drifts can lead to correlations that 
might be mistaken for intrinsic interactions. Numeri- 
cal simulations have demonstrated that the configura- 
tional temperature tracks sudden temperature jumps Q 
far faster than the pressure or internal energy, and that 
its fluctuations also decay more rapidly, particularly for 
shorter-ranged potentials [2l|. These observations sug- 
gest that the configurational temperature may converge 
to a well-defined value even before the system achieves 
equilibrium. This surprising insensitivity has been as- 
cribed [2l| to the dominant contribution to the configu- 
rational temperature of forces among nearby pairs of par- 
ticles, which can relax to near-equilibrium configurations 
long before a disturbance can propagate through the en- 
tire system. The configurational temperature, therefore, 
is a better probe of local equilibrium than global. 



D. A sum rule for interaction measurements 

Thermodynamic sum rules provide the necessary tests 
for global equilibrium. A particularly convenient form 
may be derived from Eq. Q if the pair potential is radi- 
ally symmetric. In this case, 



k B T(V 2 r u(r)) = (\V r u(r)\ 2 ) 



(29) 



and we can explicitly calculate the thermodynamic aver- 
ages of both sides of this equation: 



(V 2 u(r)) = 2%nN 



/ 1 du <i 2 u\ 

, , rg(r) dr 
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r dr dr 2 
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dr dr 



(rg(r)) dr 



f°° du dg(r) , 
2nnN / — r dr, 

dr dr 



(30) 



and 



(|V r w(r)| 2 ) = 2nnN J™ {^\ g{r)r dr. (31) 
Combining these results yields the sum rule 



du(r) dlng(r)\ du{r) 



dr 



dr 



dr 



■g(r)rdr = 0. (32) 



This sum rule should apply at arbitrary areal densities 
for any system whose interactions can be described by a 
pairwise- additive central potential, u(r). A similar result 
was obtained in Ref. (2^ for three-dimensional systems. 

Using the radial distribution function to average over 
pairs of particles removes any sensitivity to local struc- 
tural variations, and thus focuses attention on global 
properties such as the degree of equilibration. Conse- 
quently, Eq. i|32|) complements the hierarchy of configu- 
rational temperature consistency checks. 



E. Practical considerations 

1. The range of interactions 

These thermodynamic tests turn out to be exceedingly 
sensitive to imperfections in experimental data, and care 
is required to apply them meaningfully For example, 
video microscopy data necessarily is restricted to a lim- 
ited field of view, even if the sample itself is substan- 
tially larger. Particles near the edge of the field of view 
may have strongly interacting neighbors just out of sight 
whose contributions to their net force would be over- 
looked. The large apparently unbalanced forces due to 
this pernicious edge effect would grossly distort estimates 
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of the configurational temperature were they included in 
averages such as Eq. (J5J. 

To avoid this, we calculate force distributions only for 
particles whose relevant neighbors all lie within the field 
of view. Such particles lie no closer than the interaction's 
range R to the edge of the field of view. We estimate R 
from u(r) and g(r) by computing 



T(r) r |V/Mr)| 2 
= 2ir— q (r) — = -— . 

T a y ; V 2 /3u(r) ' 



(33) 



which hcuristically describes an effective contribution 
to the configurational temperature of 2irrg(r) pairs of 
spheres interacting with potential u(r) at range r [5j. A 
typical example appears in Fig. EI This should be consid- 
ered no more than a heuristic guide because it neglects 
three-particle correlations, which are known 21] to be 
important for estimating the temperature. 
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FIG. 3: Contributions to the configurational temperature due 
to pairs of colloidal spheres separated by distance r, calculated 
from the fifth-order polynomial fit to the data in Fig. Q For 
this system, spheres separated by more than R = 2<r con- 
tribute negligibly to the configurational temperature. 

Pairs at large enough separation that their interactions 
are vanishingly weak contribute negligibly to T(r). The 
interaction range R therefore can be estimated from the 
trailing edge of T(r) in Fig. By considering particles 
no closer than R to the edge of the field of view, we can 
ensure that no relevant pair interactions will be missed 
in calculating the full configurational temperature. The 
restricted field of view is indicated by the small rectangle 
overlaid on the photograph in Fig. [21 Given this cut on 
the data, we can proceed to calculate the full configura- 
tional temperature. 

The data in Fig. display two features whose validity 
we can assess using the configurational temperature for- 
malism. The first is the anomalous minimum that can 
be interpreted as evidence for long-ranged attractions be- 
tween like-charged particles. The second is the surpris- 
ingly weak contact repulsion. This low barrier to ag- 



gregation probably is not real, otherwise particles would 
aggregate rapidly by van der Waals attraction. Could 
the minimum similarly be an artifact? 

Closer inspection of the data in Fig. reveals that 
g(r) is finite even at r < a, which should be impossible. 
Two sources of experimental error, projection error due 
to particles' small out-of-plane motions and polydisper- 
sity in the spheres' diameters, artificially increase g(r) 
near contact and thus dramatically reduce the apparent 
interaction energy near r — a. What Fig. ^ shows is the 
effective potential consistent with the raw data for the 
particles' positions p(r,t), including both of these con- 
tributions. The configurational temperature calculation 
also is based on the same raw position data and so re- 
quires the associated effective potential as an input. The 
consistency condition T con fi g = 1 thus tests the accu- 
racy and thermodynamic self-consistency of u(r) for the 
measured set of p(r, t) data. 

Given p(r, t) and u(r), the net force Fj(t) on the j-th 
particle at time t can be estimated using Eq. (JSJ), with 
the sum over neighboring particles being restricted to 
those with |r< — r 3 -| < R. The set of single-particle forces 
then can be compiled into estimates for the configura- 
tional temperature for that particular snapshot, and a 
sequence of snapshots averaged to obtain a final result. 
Averaging is not necessary if each snapshot captures a 
large enough number of particles. The dilute samples 
in our study, however, typically yield N = O (100) so 
that several thousand frames are required for adequate 
statistics. 



2. Influence of measurement errors 

Even when care is taken to avoid edge effects, naively 
calculating forces, and thus temperatures, with the ex- 
perimentally sampled u(r) yields unsatisfactory results, 
as the data in Fig. 0] demonstrate. This shows the his- 
togram of T con p values obtained from the 11912 video 
frames used to generate Fig. ^ Although the distribu- 
tion is indeed peaked around unity, the average, T con p = 
1155, deviates wildly from the expected value. Other 
definitions also fare badly, with T con \ — 23, for example. 
These extraordinarily large averages are not representa- 
tive of the force distribution, however. Rather, they can 
be ascribed to a small number of frames with huge ap- 
parent temperatures. More disturbingly, many frames 
appear to have negative temperatures. 

Frames with negative temperature must have 
J2j V 2 ^ < 0. The inset to Fig. g| shows that this 
is indeed the case. Such negative values should be 
extraordinarily rare in an equilibrated system because 
they signal mechanical instability. In this case, however, 
they can be ascribed to scatter in the experimentally 
determined potential, u(r), which is greatly emphasized 
by the Laplacian operator. Particularly for r > 2a where 
u(r) is close to 0, the signal-to-noise ratio becomes small 
and errors in the computed configurational temperature 
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FIG. 4: Histogram of the configurational temperatures cal- 
culated from the measured pair potential in Fig. Q with the 
interaction range 1 < r/a < 10 for a set of 11912 video frames, 
with an average of N = 288 particles per frame. Inset: corre- 
sponding histogram of £\ V 2 iii. 




FIG. 5: Histogram of the configurational temperatures along 
the x and y directions calculated with the polynomial fit to 
the pair potential in Fig. Solid curves: T^ nF (circles) and 

TconF (squares). Dashed curve: T CO nF from raw u(r). Inset: 
Contribution to the configurational temperature due to pairs 
at separation r for the measured u(r) (dashed curve) and the 
polynomial fit to u(r) (solid curve). 



becomes unacceptably large. Cutting the long-range part 



of u(r) helps to some degree by substantially decreasing 
the number of frames with negative temperature. How- 
ever, this does not suppress the high-temperature tail in 
the histogram of single-frame temperatures, as can be 
seen from the dashed plot in Fig. [5] These artifacts still 
boost the apparent temperature to an unrepresentative 
TconF = 1150. 

To minimize the influence of scatter in u(r), we com- 
pute the configurational temperature using a fifth-order 
polynomial fit shown in Fig. ^ as the input potential. 
This simple smoothing procedure eliminates both neg- 
ative and unreasonably large values of the calculated 
configurational temperature. The solid curves plotted 
in Fig. show the resulting distribution of configura- 
tional temperatures factored along the orthogonal x and 
y directions in the field of view according to Eq. l|l(Jfl . 

The two component, T^ nF and T^/ nF differ by less than 
2.5%, which helps to confirm the system's isotropy. Their 
mean, T con p = 0.969, also is reassuringly close to the ex- 
pected value. 

Even when the temperature is calculated with the 
smoothed pair potential, a few frames still have negative 
temperatures. This is reasonable if Um r ^oo u(r) = and 
the pair potential has an extremum at an intermediate 
separation, < r < R. In this case V 2 Uj must change 
sign and some particles may contribute negative values 
to the average in a snapshot. Chances improve for the 
average itself to be negative if N is small. For example, 
10 out of 11912 frames have negative apparent tempera- 
tures in a sample with (N) = 97. Negative temperatures 
are not observed in samples characterized by purely re- 
pulsive interactions, such as the example introduced in 
Sec. EH 



3. Finite size scaling 

The relatively small number of particles in the field of 
view has other ramifications. Because the various con- 
figurational temperature definitions involve different ap- 
proximations of order 1/N, we might expect their results 
to differ from each other and from the actual thermody- 
namic temperature accordingly. Imaging a substantially 
larger region of the sample is not feasible. On the other 
hand, deliberately sub-sampling the field of view allows 
us to probe the dependence on sample size, for which we 
can extrapolate the configurational temperatures to the 
thermodynamic limit. 

The data in Fig. were obtained by calculating the 
configurational temperature with the 640 x 480 pixel 2 
field of view reduced by borders of 20, 40, 60, 80, 100, 
120, 140, 160 and 180 pixels. The interaction range in 
this system is R = la = 14.9 pixel. Reducing the number 
of particles in this manner substantially changes the esti- 
mated temperatures, thereby confirming the importance 
of finite-size scaling. We fit these results to polynomials 
in 1 /N with statistical weighting estimated from the area 
and the interaction range. These weighted fits describe 
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FIG. 6: Finite size scaling of T con F, T con \ and T con 2 for the 
data in Fig.|S] with area- weighted fits to second-, second- and 
third-order polynomials, respectively. 



the data very well, so that the extrapolation to l/N = 
should yield meaningful estimates for the configurational 
temperatures in the thermodynamic limit. Indeed, the 
extrapolated results, T con p = 0.965, T con \ — 0.956 and 
T CO n2 — 0.976, agree quite well with each other, and all 
appear to be consistent with unity. Confirming ther- 
modynamic consistency, however, requires us to assess 
the configurational temperatures' sensitivity to errors in 
u(r). 



4- Sensitivity to input potential 

The apparently good extrapolation of the configura- 
tional temperatures to the thermodynamic temperature 
would offer little insight into the nature of the confined 
colloids' interactions if these results were insensitive to 
relevant features in the pair potential. For example, 
considerable attention has been paid in the literature to 
the possibility that anomalous confinement-induced like- 
charge attractions such as the example in Fig.^might be 
artifactual. However, if we truncate the negative region 
of u(r) to create a purely repulsive potential and recalcu- 
late the configurational temperatures, T con p, T con \ and 
T con 2 all extrapolate to 1.5, an error of 150°C. The large 
deviation resulting from this admittedly crude test sug- 
gests that the observed attraction is indeed an integral 
part of the charged particles' interaction. 

More sensitive tests for thermodynamic consistency 
are provided by the hyperconfigurational temperatures 
defined in Eq. ©. Figure [7| demonstrates how small 
variations in u(r) can cause the hyperconfigurational 
temperatures to deviate with respect to each other and 
also with respect to the thermodynamic temperature. 
Two smoothed versions of the potential are plotted in 



Fig. EJa), one fit over the range 0.93 a < r < R and 
the other over the more restricted a < r < R. The 
former collapses the entire hierarchy of hyperconfigura- 



tional temperatures plotted in Fig. EJb) to ( 



1 



in the extrapolated thermodynamic limit, with Ti 
1.012, which compares favorably to T con p = 0.965 in 
Fig. EJ The latter yields the far less satisfactory results 
in Fig. [7{c). Rather than collapsing onto the thermody- 

(s) 

namic temperature, deviates systematically to lower 
values with larger index s. 

This qualitative difference is due to substantial contri- 
butions from pairs of particles with r < a. Such pairs 
should not be present in a monodisperse sample of im- 
penetrable spheres, but appear in practice because of 
the sample's 3% polydispersity in radius and because of 
projection errors due to the particles' out-of-plane fluc- 
tuations. These two effects are responsible for the ob- 
served correlations at r < er in Fig. ^ and for the un- 
reasonably small values of u(r) in the unphysical range 
0.5 a < r < a. The successful collapse of the config- 
urational and hyperconfigurational temperatures under 
these conditions demonstrates that the effective poten- 
tial accounts for the apparent particle distribution p(r) 
and may differ subtly from the ideal pair potential. 

The results in Fig. Efb) and |7fc) reflect the general 
trend that higher order hyperconfigurational tempera- 
tures are more sensitive to details of the input potential. 
Even so, we can adjust the input potential within the ex- 
perimental error bounds so that all of the hyperconfigura- 
tional temperatures converge to unity. In this sense, the 
hyperconfigurational temperatures not only strengthen 
our conclusions regarding the nature of anomalous like- 
charge attractions, but also enable us to improve our esti- 
mates for u(r) by adjusting for improved thermodynamic 
self-consistency. 

The data in Fig. [7Jc) also highlight another general 
feature of the configurational temperatures. Even though 
the more restricted trial potential does not successfully 
collapse the data, it does yield consistent results for con- 
figurational temperatures factored along orthogonal di- 
rections. This is a good indication that, indeed, the sys- 
tem is isotropic. 



5. Checking for isotropy 

Comparing temperatures measured along orthogonal 
directions not only provides a useful check for the inter- 
actions' isotropy, but also can be used to appraise the 
imaging system. For an isotropic sample, we expect that 



AT 

(T) 



(34) 



will vanish for arbitrary choices of u{r). More usually, 
the combination of a commercial video camera and video 
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FIG. 7: (a) Measured pair potential u(r) together with 
least-squares fits to fifth-order polynomials over the range 
0.93 < r/a < 2 and 1 < r/cr < 2. (b) Hyperconfigurational 
temperatures including data from r < a. = 1.012. (c) 

Hyperconfigurational temperatures over the more restricted 
range, factored into Cartesian components. 



frame grabber results in slightly different length scale cal- 
ibrations. Even quite small differences measurably affect 
the temperatures calculated from the associated force 
components if the interaction is assumed to be isotropic. 
For example, we found that AT = -0.025 (T) for the 
confined sample at H = 9 and AT = 0.06 (T) for the 
sample at H — 195 /im. Even such small differences are 
plainly visible in Fig. 0[c) . 

Apparent anisotropics of this magnitude appear con- 
sistently in our data sets regardless of the samples' com- 
position, concentration, degree of confinement, and so 
are unlikely to reflect statistical errors. Nor are they 
likely to signal a real anisotropy in our samples' interac- 
tions. Instead, they result from the hyperconfigurational 
temperatures' sensitivity to subtle geometric distortion 
in our imaging system. Rescaling the measured x and y 
coordinates slightly can substantially reduce the appar- 
ent anisotropy in the entire hierarchy of hyperconfigura- 
tional temperatures, as the data in Table[I]show. For the 
system used in this study, a 0.7% correction of the x : y 
scale ratio is enough to account for the 5.8% anisotropy 
of Xf } in the H = 195 fim data. 

The same scaling factor also corrects the apparent 
anisotropy in the other samples we have studied, and thus 
appears to be correctly interpreted as a correction to the 
calibration of our imaging system. Furthermore, differ- 
ences in the scaling factors as small as ±0.1% perform 
substantially less well, as shown by the data in Table [I] 
This level of sensitivity greatly exceeds the typical 1% 
calibration accuracy obtained by imaging test patterns, 
and thus provides a new tool for assessing and correcting 
geometric defects in the digital video microscopy system. 

Successfully correcting apparent anisotropy in the hy- 
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TABLE I: Correcting apparent anisotropy in the hypercon- 
figurational temperatures of the H = 195 /im data set by 
rescaling coordinates. In each case, g(r), u(r) and T^ s ' were 
recalculated with revised particle locations. 



perconfigurational temperatures also provides insight 
into the nature of the system's interactions. If replac- 
ing u(r) by an arbitrary function causes T x to differ from 
T y then the system indeed may be anisotropic, either 
because its interactions are anisotropic, or else because 
of its response to an external field. In the latter case, 
the external field contributes an additional configuration- 
dependent term to the Hamiltonian, Ti. ext (T), which con- 
tributes, in turn, to the definitions of the configurational 
temperatures. If d x TL ext (T) ^ d y H ex t(T), then the con- 
figurational temperatures along the two directions gen- 
erally will differ. 



6. Accounting for sample imperfections in T } 



h 



Colloidal samples quite often include small populations 
of aggregated pairs of spheres, also known as dimers. 
These can have a striking effect on the configurational 
temperature, as the data in Fig. |S] demonstrate. This 
sample also consists of a = 1.58 /im diameter silica 
spheres, but sedimented to the bottom wall of a slit pore 
H = 195 /im thick at an areal density no 2 — 0.080. 
Unlike samples in relatively thin (H — 9 /im) slit pores 
considered so far, such weakly confined silica monolayers 
are found to have purely repulsive screened-Coulomb in- 
teractions [l2iri9] | in at least qualitative agreement with 
mean-field theory. The presence of a few dimers, how- 
ever, contribute a small peak to the radial distribution 
function at r = 0.8 a whose tail extends to 1.1 a. These 
excess correlations are dramatically magnified in T(r)/T 
because any two nominally repulsive particles ought to 
exert exceptionally large forces on each other near con- 
tact. The overall result is extremely large configurational 
temperatures. 

At first glance, such defects in the sample would ap- 
pear to render meaningful measurements of the configu- 
rational temperature impossible. However, the function 
T(r)/T can be used to cut the spurious data while re- 
taining enough useful information for an accurate assess- 
ment. In particular, the peak in T{r)/T at small r is well 
separated from the principal peak at r — 1.8 a. This sug- 
gests that the former can be ascribed entirely to dimers 
and the latter to genuine long-ranged interactions, with 
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FIG. 8: Contribution to the configurational temperature at 
different pair separations assuming interactions described by 
Eq. lEU with kT 1 = 180 nm, Z = 7563 for a silica monolayer 
at H — 195 lira. The peak at 1 < r/a < 1.2 is caused by 
dimers. The tiny peak at r = 1.5 a is due to one particle 
stuck to the bottom surface. Right Inset: g(r) with a dimer 
peak around r — 0.8 a. Left Inset: a more detailed view of 
g(r) near r = a. 



a clean division at about r = 1.3 cr. It seems reasonable, 
therefore, to eliminate dimers' contribution to the config- 
urational temperatures by excluding any pair separations 
smaller than 1.3 a. 

Restricting the range of u(r) too much would exclude 
relevant interactions, causing the net force on many par- 
ticles to appear unbalanced and the configurational tem- 
perature to increase accordingly. For the present data 
set, we find that increasing the lower cutoff to r — 1.45 a 
only increases the apparent temperature by a few per- 
cent. This indicates that three-body correlations are 
weak at this concentration range and that the proposed 
cutoff at r = 1.3 cr will not distort the results. Similarly, 
ignoring contributions from pairs separated by more than 
r = 2 a has little influence on the estimated temperature 
and establishes the interaction range for accurate tem- 
perature estimates. 

Figure [5] shows typical finite-size scaling results for a 
trial fourth-order polynomial fit to the experimentally 
obtained pair potential over the range 1.435 < r/a < 
1.89. The fit potential is plotted as a dashed curve in 
Fig. HUT a). Replacing this with a fit to the predicted 
mean field potential described by Eq. (|23(l over the range 
1.2 < r/a < 1.95 yields comparably good convergence, 
and increasing the range to 1.2 < r/a < 10 decreases 
the extrapolated configurational temperature by just 4 
%. All of these trial potentials fall within error estimates 
for the measured potential. On this basis, and because 
the various definitions of the configurational temperature 
all extrapolate to unity, we conclude that the measured 
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1 /N 

FIG. 9: Finite size scaling of three configurational tempera- 
ture definitions computed for a purely repulsive monolayer of 
a = 1.58 fim silica spheres in a slit pore of height H = 195 /im. 
Solid curves are weighted fits to second-order polynomials in 
1/N used to extrapolate to the thermodynamic limit. Data at 
1/N < 0.005 are influenced by the stuck particle, as discussed 
in the text. 



potential once again accurately describes the system's 
equilibrium pair interactions. In this case, however, the 
potential appears to be purely repulsive. 

The sample used to compile these data also included 
one particle that had deposited irreversibly onto the 
lower glass surface near the edge of the field of view. 
A single immobilized sphere might not be expected to 
influence the free monolayer's structure and dynamics 
much. The influence on the radial distribution function 
is indeed subtle, with the slight peak at r = 1.5 a in 
T(r)/T (Fig.|SJl disappearing when the region containing 
the stuck particle is excluded from the calculation. The 
effect on the candidate pair potential is somewhat more 
pronounced, particularly for r < 1.5 cr, as can be seen 
in Fig. HOf a'). But does the potential from the restricted 
data set better reflect the system's interactions? After 
all, the configurational temperatures calculated with the 
unrestricted data in Fig.[5]all extrapolate reasonably well 
to the thermodynamic value. Figures llOt b) and llUf c) 
show that the hierarchy of hyperconfigurational temper- 
atures collapses to unity only when the region containing 
the stuck particle is excluded. The restricted data set, 
therefore, offers a more accurate picture of the pair po- 
tential. This is consistent with our earlier observation 
that small uncertainties in g{r) at small separations can 
be dramatically magnified in it(r). 



F. Applying the thermodynamic sum rule 

Even when the configurational and hyperconfigura- 
tional temperatures converge to the thermodynamic 
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FIG. 10: Influence of an immobilized sphere on the effective 
pair potential and hyperconfigurational temperatures, (a) 
Measured pair potentials together with fourth-order polyno- 
mial fits yielding optimal collapse of the hyperconfigurational 
temperatures. Squares: results including all data. Circles: 
results obtained by excluding the region around the stuck 
sphere, (b) Hyperconfigurational temperatures obtained for 
the entire field of view, (c) Hyperconfigurational tempera- 
tures obtained for the restricted data set excluding the region 
around the stuck particle. 




FIG. 11: Integrand of the sum rule in Eq. (1321 for the con- 
fined silica monolayer at H = 9 /mi. Inset: Best fit pair 
potential from Fig. (circles) and another estimate (squares) 
consistent with both the confidence interval of u(r) and with 
convergence of the configurational temperatures to the ther- 
modynamic value. The former gives a relative error of 0.5 in 
the sum rule, and the latter 0.001. 



value for a reasonable choice of pair potential, the sys- 
tem's degree of equilibration still has to be assessed. Fig- 
ure ^2 shows the integrand of the sum rule in Eq. I|32l) for 
two slightly different effective potentials, both of which 
are consistent with the interactions measured in the con- 
fined colloidal monolayer at H — 9 /im. Even differences 
in the input potential too small to affect the configura- 
tional temperatures can change the sum rule's integrand 
substantially. These changes affect whether or not the 
sum rule as a whole is satisfied. 

Normalizing the integral in Eq. i.'>2\) by the integral 
of the absolute value provides a useful measure of con- 
vergence. The best-fit pair potential obtained from the 
liquid structure inversion (plotted as circles in the inset 
to Fig. II ip yields an unacceptably large relative error of 
0.5. At such low areal densities (no 1 — 0.0684), u(r) 
is very similar to the potential of mean force, w(r) — 
—ksT In g(r), and very small changes in u(r) can influ- 
ence the sum rule's integrand substantially. In fact, ad- 
justing the potential just slightly to the function plotted 
as squares in Fig. llll improves the sum rule's convergence 
to 0.001. This modified potential still successfully col- 
lapses the configurational and hyperconfigurational tem- 
peratures, and thus may be considered an improved es- 
timate for the pair potential. This successful application 
of the thermodynamic sum rule suggests that the system 
is indeed in equilibrium, and that the candidate pair po- 
tential, including its long-ranged attraction, accurately 
describes the confined particles' interactions. 



G. Areal pressures 

The foregoing considerations serve to confirm that the 
same colloidal silica spheres experience qualitatively dif- 
ferent equilibrium pair interactions when confined be- 
tween parallel glass walls separated by H — 195 /im and 
H = 9 /im, with the more strongly confined dispersion 
exhibiting anomalous long-ranged attractions. Such a 
qualitative difference in the system's pair potential also 
should manifest itself in the system's other thermody- 
namic properties, such as its pressure. 

Given the particles' locations and an estimate for their 
pair interactions, we can calculate the monolayers' pres- 
sure as P = nkgT + p pot , where p 9 "* is a departure from 
ideal gas behavior due to the particles' interactions. This 
so-called "potential pressure" is give by 

\ i I i<j 

where d is the dimension of the space, = Vj — Vi 
and Fij = — V,tt(ry-). The former definition based on 
absolute coordinates works best for systems with open 
boundaries, which are described in the grand canonical 
ensemble, but fails in systems with periodic boundary 
conditions, which often are used in computer simulations. 
The latter definition based on relative coordinate works 
in both. Either should apply to our experimental data. 
Comparing the different systems' interaction-driven 
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FIG. 12: Finite-size scaling of the normalized areal pressure 
for systems with repulsive and attractive potentials. 



departures from ideal gas behavior is facilitated by defin- 
ing 



D 



(3P 



1 



fippot 



(36) 



whose dependence on system size is plotted in Fig. El 
These data were obtained with the definition based on 
absolute coordinates measured from the center of the field 
of view. Equivalent results obtained with relative coordi- 
nates differ by less than 1.5%. As for the configurational 
temperature definitions, Eq. (|35|l applies in the thermo- 
dynamic limit, and B is obtained by extrapolating to 
large system size. In this case, we see that the system at 
H = 195 /mi with long-ranged repulsive interactions has 
a substantially higher pressure than would be expected 
for an ideal gas at the same areal density. The increase is 
much smaller in the more strongly confined system, pre- 
sumably because of the potential's long-ranged attractive 
tail. 

Like the configurational temperature, the pressure also 
depends sensitively on the input potential and length 
scale calibration. Without appropriate rescaling, P x is 
10% higher than P y in the H = 195 /im data set. After 
rescaling using the factors obtained by requiring isotropy 
in the configurational temperatures, the two components 
agree well, as can be seen in Fig. This further con- 
firms that the observed anisotropy is due to artificial im- 
age distortion that can be corrected with a single rescal- 
ing factor. 



H. 



Using to measure u(r) 



The configurational and hyperconfigurational tem- 
peratures show great promise as thermodynamic self- 
consistency tests for interaction measurements in soft- 



matter system. They also can be used to measure interac- 
tions directly, with the hierarchy of definitions providing 
a set of simultaneous constraints on free parameters in a 
model for the pair potential. Even when the theoretical 
form of the potential is not known a priori, a numeri- 
cally determined function that simultaneously collapses 
the entire hierarchy to the thermodynamic temperature 
would provide a model-free estimate for the potential. 

As an example of this procedure, we use hypercon- 
figurational temperatures to determine the parameters 
k _1 and Z in the screened- Coulomb potential, Eq. H23|) . 
that accounts for the distribution of particles in the 
H = 195 /im data set 0. Requiring = T^ 1 = 1 
yields a screening length of k _1 = 208.4 nm and an ef- 
fective charge on the particles of Z = 4123. These values 
also converge and close to unity, suggesting that 
the screened- Coulomb potential reasonably describes the 
interactions in this system. Table ITU shows typical results 
for parameters near the optimal values. 

Simply fitting Eq. (|23|l to the measured pair po- 
tential in Fig. HUT al yields k _1 = 180 ± 10 nm and 
Z = 6500 ± 1000. Although these values differ some- 
what from those obtained by collapsing the configura- 
tional temperatures, both are broadly consistent with the 
charge regulation theory for interacting silica surfaces |9j . 
Their disagreement and the slight departure from unity 
of the highest-order hyperconfigurational temperatures 
suggests a small deviation of the true pair potential from 
the predicted screened-Coulomb form. This is reasonable 
because the nearby bottom wall probably affects the in- 
teraction somewhat. No convergence is possible at all 
for the data at H = 9 /im because the pair interactions 
are not well described by the purely repulsive screened- 
Coulomb form. 

Using hyperconfigurational temperatures to provide 
constraints on models for the pair potential offers ad- 
vantages over other interaction measurement techniques. 
Principally, this method eliminates the intermediate 
steps of first measuring g{r) and then inverting it to 
obtain u(r). This first step generally requires amassing 
large amounts of data to obtain accurate information on 
correlations at small separations. The second involves 
uncontrolled approximations and limits the range of con- 
centrations over which interactions can be measured. By 
contrast, the present approach requires substantially less 
data and involves fewer approximations. Consequently, 
this approach should be useful for exploring denser, more 
strongly interacting systems. Its ability to distinguish 
pairwisc and many-body interactions also should be help- 
ful for such studies. 

On the other hand, calculating configurational temper- 
atures involves sums over all pairs of particles, and must 
be repeated for sub-samples of varying size to account 
for finite-size scaling. This can be computationally ex- 
pensive, particularly in model-free searches. Implement- 
ing fast search algorithms on restricted data sets with 
additional data being added only in the final stages of 
polishing should alleviate this problem. 
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K 1 (nm) 


Z 


1 h 


itt(3) 

M h 


m(5) 

1 h 


rp(~) 

1 h 


160 


10147 


1.00 ±0.03 


1.62 ±0.12 


1.53 ± 0.09 


1.23 ±0.06 


180 


6755 


1.00 ±0.02 


1.31 ±0.09 


1.31 ±0.07 


1.08 ±0.05 


200 


4779 


1.00 ±0.01 


1.07 ±0.07 


1.12 ±0.06 


0.94 ± 0.04 


210 


4096 


1.00 ±0.01 


0.98 ± 0.06 


1.04 ±0.06 


0.88 ± 0.04 


220 


3549 


1.00 ±0.01 


0.90 ±0.05 


0.96 ± 0.05 


0.82 ± 0.04 


250 


2438 


1.00 ±0.01 


0.72 ± 0.03 


0.77 ±0.04 


0.68 ±0.03 



TABLE II: Determining the free parameters in a trial potential using . Hyperconfigurational temperatures calculated from 
the screened- Coulomb potential for the H — 195 /im data set over the range 1.3 < r/cr < 10. 



Failure of such a search to converge may signal a failure 
of pairwise additivity, a departure from equilibrium, or 
else an inappropriate model for the pair potential. Dis- 
tinguishing these would probably require resorting to one 
of the alternate measurement methods. 



VI. CONCLUSION AND DISCUSSION 

The recently introduced notion of configurational tem- 
perature provides a valuable new tool for assessing exper- 
imental systems' thermodynamic state from static snap- 
shots. We have introduced a hierarchy of hyperconfigura- 
tional temperatures that emerge naturally from the gen- 
eralized definition of temperature, and have shown that 
both these and the generalized definition of the temper- 
ature can be derived from the classical hypervirial theo- 
rem. Colloidal monolayers provide an ideal experimental 
test bed for these new concepts in statistical mechanics. 

We have tested the configurational and hypercon- 
figurational temperature to within 1% accuracy using 
particles' distribution and pair potentials measured in 
Ref. 0|- The effect of finite system size is clearly 
observed and accounted for in these measurements. 
Since the configurational temperatures' derivation re- 
quires pairwise additivity and their computation depends 
sensitively on the input potential, they can be used as 
thermodynamic self-consistency checks for measured pair 
potentials. The configurational temperatures calculated 
for our experimental data on confined colloidal silica 
monolayers confirm that our measured potentials are ac- 
curate and that they are consistent with the assumption 
of pairwise additivity. 

We have found that higher-order hyperconfigurational 
temperatures are increasingly sensitive to errors in the 
potential u{r) because of their dependence on higher mo- 
ments of the force distribution. Even so, we have found 
that we can adjust the pair potentials within the mea- 
sured uncertainties to converge the entire hierarchy of 
hyperconfigurational temperatures to unity for all of our 
data sets. This provides substantial independent evi- 
dence that our directly measured potentials reliably re- 
flect equilibrium pair potentials for our systems. 

A sum rule introduced in Sec. IV Dl complements the 
information provided by the configurational and hyper- 



configurational temperatures by providing additional in- 
sights into the system's degree of equilibration. Very 
small errors in the pair potential, moreover, are dramati- 
cally emphasized by the derivatives in the sum rule's def- 
inition. We find, nevertheless, that the sum rule can be 
satisfied for the colloidal samples in our study by adjust- 
ing the potential within the range of their uncertainties. 

Applying these analytical tools to our system of sedi- 
mented colloidal silica spheres allows us to draw new con- 
clusions regarding the nature of electrostatic interactions 
in this deceptively simple system. Rather than being 
purely repulsive, as mean-field theory predicts, confined 
colloids' interactions can be characterized by a strong and 
long-ranged attraction. This result echos those reported 
more than a decade ago in the first generation of colloidal 
interaction measurements [HQ 111 HI HI Now, how- 
ever, we can assert with confidence that the observed 
anomalous attractions constitute pairwise-additive con- 
tributions to the systems' equilibrium free energies, and 
not from any of the myriad of possible experimental 
artifacts that have been proposed. This interpretation 
is further bolstered by measurements of the areal pres- 
sure in these monolayers, which show clear signatures 
of their differing interactions. Resolving the anomaly of 
confinement-induced like-charge attractions therefore re- 
quires a fresh assessment of the nature of colloidal elec- 
trostatic interactions in simple electrolytes. 

Hyperconfigurational temperatures also can be used as 
a set of constraints to determine the free parameters in a 
model for a system's pair interactions. We have demon- 
strated this by determining the two free parameters in 
a screened- Coulomb model for charge-stabilized colloids' 
interactions in the repulsive regime. In principle, the un- 
bounded hierarchy of hyperconfigurational temperatures 
can be used in this way to determine any pairwise ad- 
ditive potential given the position, particularly if that 
potential can be modeled as a polynomial or comparably 
simple function. 

One advantage of this method is that it circumvents 
the uncontrolled approximations that have dogged other 
approaches to measuring macroionic interactions in equi- 
librium. This method in more general in that it can 
be applied at arbitrary particle densities. Unlike mea- 
surement techniques based on the radial distribution 
function, g(r), furthermore, configurational temperature 
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measurements can be applied to inhomogeneous systems. 
This new method therefore should be useful in phase sep- 
arated systems at equilibrium. 

We are grateful to Owen Jcpps, Sven Behrens and 
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